Submitted to JGR - Atmospheres Aug 2000 


Northern Hemisphere winter climate response to greenhouse gas, 
ozone, solar and volcanic forcing 


Drew T. Shindell 

NASA Goddard Institute for Space Studies, New York 
Gavin A. Schmidt 

NASA Goddard Institute for Space Studies and Center for Climate Systems Research, 
Columbia University, New York 


Ron L. Miller 

NASA Goddard Institute for Space Studies, New York 
David Rind 

NASA Goddard Institute for Space Studies, New York 


Abstract 

The Goddard Institute for Space Studies (GISS) climate/middle atmosphere model has been 
used to study the impacts of increasing greenhouse gases, polar ozone depletion, volcanic 
eruptions, and solar cycle variability. We focus on the projection of the induced responses onto 
Northern Hemisphere winter surface climate. Changes in the model’s surface climate take place 
largely through enhancement of existing variability patterns, with greenhouse gases, polar ozone 
depletion and volcanic eruptions primarily affecting the Arctic Oscillation (AO) pattern. 
Perturbations descend from the stratosphere to the surface in the model by altering the 
propagation of planetary waves coming up from the surface, in accord with observational 
evidence. Models lacking realistic stratospheric dynamics fail to capture these wave flux changes. 
The results support the conclusion that the stratosphere plays a crucial role in recent AO trends. 
We show that in our climate model, while ozone depletion has a significant effect, greenhouse gas 
forcing is the only one capable of causing the large, sustained increase in the AO observed over 
recent decades. This suggests that the AO trend, and a concurrent strengthening of the 
stratospheric vortex over the Arctic, are very likely anthropogenic in origin 



2 


1. Introduction 

Recent observations show very large increases in 
surface temperatures over the Northern Hemisphere 
continents during winter (Figure lA). These result 
largely from an increase in the westerly air flow 
around the Northern Hemisphere, which brings warmer, 
wetter oceanic air over the continents. This westerly 
flow, and therefore the surface temperature changes, 
are associated with the leading variability pattern of 
Northern Hemisphere cold-season sea-level pressure 
(SLP), called the ’Arctic Oscillation’ (AO) (Figure 
IB). Observations show an apparent upward trend 
in the amplitude of this pattern, or equivalently, a 
bias towards the positive phase of the pattern, in re- 
cent decades [Thompson and Wallace, 1998; Thomp- 
son et ai, 2000]. The variability associated with this 
pattern extends from the surface up into the strato- 
sphere, and in fact the variability pattern can be 
equivalently defined as being composed of variability 
at all levels from the surface through the lower strato- 
sphere [Baldwin and Dunkerton, 1999]. Previously, 
we showed that only with the inclusion of a well- 
resolved stratosphere in our climate model were we 
able to reproduce the observed trend in the AO in re- 
sponse to increasing greenhouse gases [Shindell et ai, 
1999a]. Consistent with those results, observations in- 
dicate that changes in the stratosphere typically pre- 
cede the changes in the lower atmosphere [Kodera and 
Koide, 1997; Baldwin and Dunkerton, 1999]. As a re- 
sult, the stratospheric climate model was better able 
to reproduce the large continental wintertime warm- 
ing trends seen in the Northern Hemisphere than a 
similar ’tropospheric model’. 

AO trends in the earlier climate model simula- 
tions were induced by the effect of increasing green- 
house gases on stratospheric temperatures. Other 
perturbations to the stratosphere also affect surface 
climate. For example, it has been shown that the 
circulation in the lower atmosphere can be affected 
by stratospheric changes such as large volcanic erup- 
tions [Rind et al, 1992; Graf et ai, 1993; Kodera, 
1994] or solar forcing [Haigh, 1996, 1999; Shindell 
et al, 1999b]. We now investigate in more detail the 
mechanism by which stratospheric perturbations af- 
fect surface changes. Greenhouse gas increases, vol- 
canic aerosols, ozone depletion or solar forcing can 
all induce meridional temperature gradients in the 
stratosphere during wintertime. We demonstrate a 
positive feedback whereby zonal wind anomalies in- 
duced by these thermal gradients near the tropopause 
are amplified through planetary wave refraction. De- 


flection of upward-propagating tropospheric waves at 
lower and lower altitudes in effect carries the anomaly 
steadily down from the stratosphere, allowing strato- 
spheric changes to affect surface climate by altering 
tropospheric energy flow. 

Additionally, we examine the extent to which each 
individual stratospheric forcing excites with the nat- 
ural patterns of variability. This is accomplished by 
projecting the induced changes onto empirical orthog- 
onal functions (EOFs), fixed spatial patterns of vari- 
ability ranked by the amount of variability accounted 
for by each pattern [e.g. Kutzbach, 1970], derived from 
a control run with no external forcings. Examina- 
tion of the observations reveals that the recent trend 
in sea-level pressure has occurred primarily through 
enhancement of the AO (EOF 1) [Thompson et al, 
2000] and Figure 2, discussed further below). In or- 
der to reproduce the observations, it’s therefore nec- 
essary that a model yield a trend in the amplitude 
of the AO pattern similar to the observations, while 
at the same time keeping the amplitudes of changes 
in the other variability patterns small. Analysis of 
our model simulations shows that while greenhouse 
gas increases, volcanic aerosols, or Arctic ozone de- 
pletion can excite primarily this leading mode, only 
increasing greenhouse gases can excite an AO trend 
comparable to the observed value 

2. Experimental Setup 

Simulations of the climate response to increasing 
greenhouse gases, Arctic ozone depletion, volcanic 
eruptions, and solar cycle variability have been per- 
formed (Table 1). These are the primary forcings 
that affect the stratosphere, which seems to play an 
important role in Northern Hemisphere winter cli- 
mate change [Shindell et ai, 1999a]. Increases have 
been observed in the abundance of atmospheric green- 
house gases, Arctic winter/spring ozone loss, and so- 
lar irradiance over decadal and longer time scales [e.g. 
Houghton et o/., 1992]. Solar cycle variability is used 
here as a proxy for longer-term solar variations, which 
are not well known. While there is not thought to 
have been a long term increase in volcanic eruptions, 
those simulations complete the range of forcings, and 
therefore allow us to perform a systematic compari- 
son of the mechanism by which perturbations to the 
stratosphere affect surface climate. 

The primary model used here is an older version 
of the Goddard Institute for Space Studies (GISS) 
stratospheric climate model which has 8° x 10° de- 
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gree horizontal resolution, and 23 vertical layers ex- 
tending up into the mesosphere (top at 0.002 mb, 
km) [Rind et al, 1988b, a]. While a newer version of 
the model, with finer resolution, is available, compu- 
tational resources do not yet permit us to carry out 
simulations with that model for the large number of 
years used here. Note that the propagation of plan- 
etary waves is a largely quasi-geostrophic process, so 
it should not be greatly affected by resolution, espe- 
cially for the longest waves which are resolved by even 
coarse resolutions. However, the amplitude of the 
waves may change with resolution, leading to quanti- 
tative differences although the results should remain 
qualitatively the same. 

This atmospheric model is coupled to a mixed-layer 
ocean, allowing sea-surface temperatures to respond 
to external forcings. A gravity wave parameterization 
is employed in the model in which the temperature 
and wind fields in each grid box are used to calculate 
gravity wave effects due to wind shear, convection, 
and topography. Planetary scale waves are generated 
by land-ocean contrasts and by the influence of topog- 
raphy. For the runs described here, the mountain drag 
was set to one quarter the amount used in the origi- 
nal documentation [Rind et a/., 1988b, a], a standard 
modification used in many of our subsequent stud- 
ies, especially those relating to polar ozone [Shindell 
et ai, 1997, 1998a, b, 1999a]. The altered model now 
gives a good reproduction of observed temperatures 
in the polar lower stratosphere of each hemisphere 
during their respective winters and springs. As in the 
original control simulations [Rind et al, 1988b, a], 
the residual circulation and the generation of plan- 
etary waves in the altered model are in fairly good 
agreement with observations. Reducing the drag does 
marginally reduce the variability in the model, as well 
as slightly degrading the simulation at other seasons 
and latitudes, however. 

The model was run with increasing greenhouse 
gases starting in 1959, about the earliest time we have 
good quality global data with which to initialize the 
model, and extending for 110 years. Given that green- 
house gas increases began much earlier, the model will 
lack the unrealized warming caused by those earlier 
emissions (the ”cold start” problem), and will there- 
fore slightly underestimate the impact of increasing 
greenhouse gases. Greenhouse gas trends were based 
on observations through 1985, and subsequently fol- 
lowed IPCC projections [Houghton et al, 1992]. Five 
versions of this run were completed, a control run 
without greenhouse gas increases, and four exper- 


imental runs. Three of the experimental runs in- 
cluded additional changes: one included changes in 
polar lower stratospheric ozone, as in Shindell et aL 
[1998a, 1999a], while the second and third both had 
changes in stratospheric ozone outside the polar re- 
gions, with the third also including positive trends 
in stratospheric water vapor amounts, which affects 
both temperatures and stratospheric ozone chemistry 
[Shindell, 2000]. 

To test the impact of volcanic aerosol perturba- 
tions on climate variability patterns, we ran the model 
again from 1959 to 2000 with increasing greenhouse 
gases, and also including stratospheric aerosol load- 
ing based on observations [Sato et al, 1993, updated 
through 1999]. The opacity of the volcanic aerosols 
varied as a function of time; winters with an opacity 
greater than .02 were deemed volcanic years (1964- 
66, 1969-70, 1975, 1983-85, 1991-94, mean opac- 
ity = .040), while years with a lower opacity were 
considered background (non-volcanic) years (1959- 
60, 1967-68, 1971-74, 1976-82, 1986-90, 1995-2000, 
mean opacity = .009). Both the volcanic and back- 
ground years are centered around 1979, so that the 
difference between the average response of the two 
sets should reflect only the volcanic forcing, and not 
the influence of the changing greenhouse gases over 
time. Secondary interactions between ozone and vol- 
canic aerosols are not included in this simulation, 
though the aerosols do influence local temperatures. 
Note that the peak opacity after major eruptions can 
be quite large, with large variations from eruption to 
eruption: for example, after El Chichon it was about 
15 times larger than the mean opacity of the non- 
volcanic years, and after Mt. Pinatubo it was roughly 
25 times larger. Use of many volcanic years allows us 
to improve the statistical power of our analysis. 

Severe ozone depletion in the polar regions, ’ozone 
holes’, were introduced into the climate model using 
an interactive parameterization of the relevant chem- 
istry, as described previously [Shindell et al, 1997]. In 
brief, chlorine activation (conversion from reservoir to 
reactive molecules) is calculated in each model grid 
box. Complete activation takes place in five hours, 
based on a simple activation switch triggered when- 
ever temperatures fall below 195°K, a rough thresh- 
old for polar stratospheric cloud (PSC) or supercooled 
ternary solution (STS) formation. This simple acti- 
vation switch gives reasonable results since the het- 
erogeneous reactions that activate chlorine are very 
rapid even for very small surface areas. When hetero- 
geneous processing ceases, deactivation follows two- 
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dimensional model derived rates. Ozone depletion 
takes place at each point in the GCM where there 
is both active chlorine and sunlight, which is required 
for the chlorine catalytic cycle of ozone destruction. 
We also include an additional estimated contribution 
to halogen catalyzed ozone depletion of 15% from 
bromine, as calculated in the offline chemical model. 
Ozone recovery rates are parameterized so that ozone 
losses are restored based on the photochemical life- 
time of ozone at each altitude. This parameteriza- 
tion hcis been shown to yield a good reproduction 
of observed levels of ozone depletion in both hemi- 
spheres during the 1980s and 1990s [Shindell et a/., 
1997, 1998a, b]. This model was run for 19 years 
with no changes in external forcing, and total chlo- 
rine set to 3.2 parts per billion by volume (ppbv), 
corresponding to early 1990s abundances. 

While long-term changes in solar irradiance have 
certainly influenced the stratosphere, they are not 
well constrained due to the lack of observations. The 
roughly 11-year solar cycle has been well observed 
however, so we have investigated its impact in our 
model [e.g. Balachandran and Rind^ 1995; Rind and 
Balachandran, 1995; Shindell et al., 1999b]. The at- 
mosphere’s response to longer-term solar output vari- 
ability would perhaps be similar, though this could 
be affected by the longer time scale of the forcing 
or by differences in the spectrum of irradiance vari- 
ability. One should keep in mind that there is also 
no guarantee that non-linear changes in atmospheric 
circulation would extrapolate linearly to longer-term 
solar changes. The simulations analyzed here in- 
clude realistic wavelength-dependent solar cycle ir- 
radiance changes from 180 to 400 nm, and constant 
changes at longer wavelengths consistent with total 
solar cycle irradiance variations [Lean et al, 1997]. 
We also includes an interactively calculated ozone re- 
sponse to those changes. The ozone response is based 
on a two-dimensional model-derived chemistry pa- 
rameterization which includes wavelength-dependent 
ozone response to changes in radiation and response 
to changes in temperature [Shindell et ai, 1999b]. So- 
lar variability directly affects both ozone photochem- 
istry and local heating, modifying ozone abundances, 
which in turn further alter local heating rates, as well 
as the radiation field at other levels. The GCM was 
run for 20 years each at constant solar maximum and 
solar minimum conditions to get reasonable statistics 
on each phase of the solar cycle. Since either phase 
would not actually persist continually for so long, sea- 
surface temperatures were held fixed in these runs so 


as not to yield a false response. 

Most calculations of the climate response to an- 
thropogenic greenhouse gas forcing have been based 
upon tropospheric models with limited stratospheric 
resolution, and a highly damped or truncated strato- 
spheric circulation. For comparison with the strato- 
spheric model, we have examined simulations with a 
tropospheric model having the same 8° x 10° horizon- 
tal resolution and physics (model II) [//arisen et a/., 
1983], and a newer version with 4° x 5° horizontal res- 
olution and several improvements in model physics 
(model IT) as discussed or referenced in Rind and 
Lemer [1996]. Both versions have 9 vertical layers, 
with only 1 or 2 layers in the stratosphere (depend- 
ing upon latitude). A simulation with a version of 
model IF containing an interactive simulation of sul- 
fate chemistry has been analyzed as well [Koch et al, 
1999]. Lastly, we have also examined the behavior 
of a coupled ocean-atmosphere version of the model 
[Russell et ai, 2000]. Both versions of the atmo- 
spheric GCM, as well as the coupled model, have been 
run with increased greenhouse gas forcing as for the 
stratospheric model. For sulfate, we examine the cli- 
mate response to anthropogenic sulfate increases by 
looking at the difference between a run with present 
day sulfate and one with pre-industrial values [Koch 
et a/., 1999]. 

We point out that 650 years in total were simulated 
with the stratospheric model, reemphasizing why we 
were forced to use a coarse resolution model 

3. Analysis 

Empirical Orthogonal Functions (EOFs) were cal- 
culated from the NH cold season (Nov- Apr) SLP time 
series of the control run. The trend from the tran- 
sient greenhouse gas runs was then projected upon 
those EOFs. (However, since the spatial pattern of 
the leading mode of SLP variability is a robust fea- 
ture resembling the AO in each simulation, we would 
have obtained similar results using each simulations 
own leading mode, as in Shindell et al. [1999a]). Other 
forcings were not run as transients, so differences be- 
tween two states are used. For the solar cycle vari- 
ability simulations, the difference between the runs 
with solar maximum and solar minimum conditions 
was used. For the volcanic forcing, the difference be- 
tween the years when stratospheric aerosol loading 
was large - primarily the few years immediately fol- 
lowing eruptions - and years with background levels 
was used. Results for the ozone hole simulations are 
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differences between the simulation with polar hetero- 
geneous chemistry included and a control run with- 
out the chemistry parameterization. For the transient 
simulations, we define the time series corresponding 
to the leading EOF as the ’index’ of the model’s AO. 
The spatial pattern has been normalized so that the 
value of the AO index is equal to the SLP anomaly 
averaged poleward of 60°N, and is given in mbar. For 
the transient runs, the AO index trend is given over a 
thirty year period to match the observational record. 
For the stratospheric model, this is the period after 
initial spin up, which requires roughly 20 years for 
the ocean to fully respond to atmospheric forcing, and 
before an apparent saturation (see discussion below). 
While this does not match exactly in time with the ob- 
servational period, as our models runs began in 1959, 
so spin-up was completed only in 1979, it seems to be 
the most realistic period for comparison, as the real 
world does not have an initial ’cold start’ problem. If 
we were to compare the same 1967-1997 period with 
the observations, the model’s trend would be reduced 
by about 1/3, but is still much larger than that seen 
in the tropospheric models (as is true for the entire 
run length, as shown in Shindell et al. [1999a]), The 
three- decade value we present here Wcis consistent for 
the entire five decades during which the AO increased 
steadily between the initial spinup and eventual sat- 
uration. For the tropospheric runs, the averaging pe- 
riod has a minimal effect on the results, as the trends 
were so small that no saturation was seen. For the 
solar, volcanic, and ozone hole simulations, whose re- 
sults are differences, we obtain a single value for the 
change in the AO index based upon the component 
of the total SLP difference associated with the AO. 

While the leading EOF is well separated from sub- 
sequent patterns, the higher EOFs are not adequately 
distinct from one another. In the stratospheric model 
forced with increasing greenhouse gases, EOF 1 ac- 
counts for 13.3% of the variance, but subsequent 
EOFs account for 5.9, 4.8, 4.1, and 3.5%, respectively, 
which are not statistically different within the calcu- 
lated variance range of about 2% (calculations based 
on the entire 110 year run). An EOF analysis of the 
observations yields similar results, with 21.3, 11.7, 
11.1, 8.2, and 7.0% variance explained by the first 
five EOFs (calculations based on 1947-1997 data). 
Again, only the leading pattern is well-separated from 
the others. Since the higher modes are mathematical 
constructs which are forced to be orthogonal to the 
leading mode, not well separated, and lacking an ob- 
vious physical interpretation, we restrict our analysis 


to the leading EOF (the AO). 

As discussed in Shindell et al. [1999a], the version 
of the GISS model containing a realistic representa- 
tion of stratospheric processes produces an increasing 
AO trend in response to greenhouse gas increases of 
?^0.8 mb/decade during the several decades between 
spin-up and saturation during which the AO index is 
increasing steadily. This value is comparable to the 
roughly 1 mb/decade suggested by observations in re- 
cent decades. In contrast, model versions lacking a re- 
alistic stratosphere, ’tropospheric models’ with only 
1 or 2 layers in the stratosphere, produce a leading 
EOF that also looks very much like the AO, but its 
index shows only a weak increasing trend (+0. 1-0.3 
mb/decade) which is not statistically significant (Fig- 
ure 3). It is interesting that when stratospheric wa- 
ter vapor is allowed to increase steadily, as a result 
of tropopause warming and methane oxidation [Shin- 
dell, 2000], this distinction between stratospheric and 
tropospheric models is slightly reduced. When water 
vapor and greenhouse gases both increase, the AO 
trend is reduced to «0.6 mb/decade, outside the range 
of the trend in the other three increasing greenhouse 
gas simulations (0.83, 0.78, and 0.79). This is due 
to the slight decrease in the meridional temperature 
gradient that results from mid-latitude cooling that 
ensues from water vapor induced ozone destruction in 
the lower stratosphere [Shindell, 2000]. Those simu- 
lations do not include the potential effect of increased 
water vapor on heterogeneous ozone destruction in 
the polar region, however, which would oppose the 
mid-latitude ozone induced cooling and could poten- 
tially cause a significant increase in the meridional 
temperature gradient in the spring (see section 5.2 
below). Water vapor on its own, as a greenhouse gas, 
increases the bias towards the high phase of the AO. 

We find that the inclusion of an ozone hole in the 
model, another candidate for driving the observed 
increases [e.g. Graf et al., 1998; Volodin and Galin, 
1999], does cause an increase in the model’s AO index, 
but of only «0.6 mb (1990s ozone depletion relative to 
the 1970s and earlier, when there was no depletion, 
or approximately 0.3 mb /decade over the past two 
decades). Though this has likely contributed to the 
measured trend, it seems insufficient to account for 
a great deal of the observed increase of several mbar 
over the past three decades. This is perhaps not sur- 
prising, since the trend in the AO has been observed 
throughout the November-to-April cold season, but 
severe ozone depletion in the Arctic does not take 
place in general until the spring, when sunlight falls 
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on cold vortex air. These results are consistent with 
our earlier conclusion that the addition of ozone de- 
pletion to greenhouse gas forcing was not required to 
induce a trend in the AO [Shindell et a/., 1999a]. This 
conclusion is also in agreement with the mild increase 
in the wintertime meridional temperature gradient at 
Northern high latitudes calculated based upon TOMS 
ozone trends [e.g. Ramaswamy et al.^ 1996]. While 
the TOMS observations were extrapolated to high 
latitudes during the polar night, the use of in situ 
ozonesonde data [Randel and Wu^ 1999a] does not 
change this conclusion [Rosier and Shine ^ 2000]. This 
contrasts with the situation in the Antarctic, where 
ozone depletion seems to account for nearly the en- 
tirety of the observed temperature trends [Randel and 
Wu, 1999b]. 

Volcanic forcing can be strong in the years imme- 
diately following an eruption, leading to an AO dif- 
ference of 0.4 mb between volcanic and non-volcanic 
years, but since it is an intermittent forcing that de- 
cays rapidly, it also seems unlikely to have contributed 
greatly to the long-term observed trend. Solar cycle 
variability causes an increase in the AO index of 0.3 
mb between solar maximum and minimum. Since the 
longer-term trend in irradiance over the past 30 years 
has been of comparable size to solar cycle variability, 
it also seems unlikely that solar variability has been 
responsible for much of the observed trend. Further- 
more, estimated solar irradiance increased as much 
in the first half of the 20th century as in the second 
[Lean et a/., 1997], but the AO index showed no in- 
crease during the former. Note that solar irradiance 
over the satellite period is much better know than 
earlier. 

The stratosphere model is distinct from the tropo- 
spheric models not only for exhibiting a trend of mag- 
nitude comparable to the observed value, but because 
the spatial distribution of its SLP trend closely re- 
sembles the AO spatial pattern, as observed [Thomp- 
son and Wallace, 1998]. The resemblance of the SLP 
trend to the AO spatial pattern can be quantified by 
decomposing the trend into the EOFs of SLP, and 
measuring the contribution of each EOF to the trend 
variance (Figure 2), following Fyfe et al. [1999]. The 
trend in observations of cold season SLP between 1947 
and 1997 [TVenberth and Paolino, 1980, updated to 
1997] is dominated by a change in the AO (the lead- 
ing EOF), which contributes 56% of the trend. The 
GISS stratospheric model with increasing greenhouse 
gases puts 64% of the trend in the AO (the average 
value of the ensemble of four simulations). The strato- 


spheric model with either an ozone hole or volcanic 
forcing also projects primarily onto the first EOF, 
however solar cycle variability does not. The strato- 
spheric model with increasing water vapor projects 
less strongly onto the AO than do the other increasing 
greenhouse gas runs, but its trend is still dominated 
by the AO (and actually looks most like the observed 
trend projections). 

In contrast, GISS model versions lacking a detailed 
stratosphere - ’tropospheric models’ - forced with in- 
creasing greenhouse gases produce trends that project 
only 5-20% onto the AO (Figure 2). The GISS tro- 
pospheric model including sulfate aerosols puts just 
6% of its trend into the leading EOF, with most of 
its trend spread out over a large number of patterns. 
Transient greenhouse gas experiments by other mod- 
eling groups seem to be consistent with this distinc- 
tion (Figure 2). The Canadian Climate Center CC- 
Cma model [Fyfe et ai, 1999], which shows a positive 
AO trend in response to increasing greenhouse gases 
despite the absence of detailed stratospheric dynam- 
ics, similarly fails to reproduce the observed domi- 
nance of the AO pattern, putting less than 30% of 
its total trend into its leading EOF. A Geophysical 
Fluid Dynamics Laboratory (GFDL) climate model 
version lacking a detailed stratosphere also finds that 
the majority of its SLP trend in response to increas- 
ing greenhouse gases is not in the AO (Paul Kush- 
ner, personal communication). Several other groups 
have presented the simulated AO response to increats- 
ing greenhouse gases. The ECHAM3 model from the 
Max Planck Institute for Meteorology, without a full 
representation of the stratosphere, shows a system- 
atic increase in the AO index in response to increas- 
ing greenhouse gases [Graf et al, 1995, 1998; Perlwitz 
et al, 2000], but its amplitude is stated to be weaker 
than observed. Paeth et al [1999] present results of 
simulations with both the ECHAM3 and ECHAM4 
models (showing NAO trends, but these are likely 
very similar to AO trends). Zorita and Gonzalez- 
Rouco [2000] present the AO trends from the same 
simulations, adding in results from the Hadley Cen- 
ter model 2 as well. Along with Fyfe et al [1999], all 
of these groups use normalized units to present their 
trends, making it impossible to compare the magni- 
tude of the increase in those simulations to the ob- 
served value. We suggest that documentation of a 
model’s AO trend should include, at minimum: (1) 
the trend’s magnitude in mb and its statistical signif- 
icance, with a description of the area weighting used, 
and (2) the decomposition of the trend into its com- 
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ponent EOFs. This information would facilitate a 
comparison across models to document more clearly 
the role of stratospheric dynamics. As yet, only mod- 
els that realistically simulate the stratosphere have 
quantitatively demonstrated an AO trend of roughly 
the right magnitude and with the same predominance 
within the SLP trend as the observations 

4. Impacts 

Strong trends in surface winds over the North At- 
lantic are closely correlated with the AO, which has 
the largest effect in this sector. Figure 4 shows trends 
in zonal wind speed from the NCEP reanalysis data 
set [Kalnay et a/., 1996] and from various GISS mod- 
els. Clearly the stratospheric model is best able to 
simulate the observed increase in westerly flow (note 
that like many climate change signals, this one is just 
beginning to emerge from an extremely variable pat- 
tern). The consequences are visible downstream over 
Siberia, where an increase in warm oceanic air carried 
from the Atlantic leads to a rise in surface air temper- 
atures. Again, the stratospheric model is best able to 
reproduce the observed changes (Figure 5 - see also 
Shindell et al [1999a]). 

Another region of great interest is the area around 
Labrador and Newfoundland, where surface tempera- 
tures have been decreasing over recent decades. Sur- 
face wind changes in this region include a significant 
meridional component, which is associated with al- 
tered planetary wave propagation, in addition to the 
strengthened westerlies discussed above. In the tro- 
pospheric model, this region warms significantly as 
greenhouse gases increase (Figure 6). The strato- 
spheric model, with an increased AO bringing more 
cool, continental air to the Eastern edge of North 
America, does a better job, with much less warm- 
ing, but still does not capture the observed cooling 
trend. In this region, the use of a coupled, dynamic 
ocean makes a considerable difference, however, as 
shown in Figure 6. The tropospheric model coupled 
to a dynamic ocean model [Russell et al., 2000] does 
a much better job than the tropospheric model run 
with a simple mixed-layer ocean. This may indicate 
that both the stratosphere and the ocean play im- 
portant roles in bringing about the cooling seen in 
this region. Additionally, the tropospheric model, us- 
ing fixed sea-surface temperatures, also does a bet- 
ter job when sulfate aerosol changes are incorporated, 
though the ones included here are pre-industrial ver- 
sus present-day, and so much larger than the changes 


over the past few decades alone 

5. Physical Mechanism 

5.1 Effect of increasing greenhouse gases 

The radiative impact of greenhouse gases at a par- 
ticular location depends upon the ratio of the absorp- 
tion of upwelling longwave radiation and the emission 
of radiation by the molecules. Since absorption and 
emission depend upon the fourth power of temper- 
ature, the net balance will be determined by differ- 
ences between the mean temperature at which the 
upwelling radiation was emitted and the local tem- 
perature at which the greenhouse gases are emitting. 
This leads to opposing radiative effects in the tro- 
posphere and stratosphere. In the troposphere, the 
lapse rate is positive (temperature decreasing with 
altitude), so that molecules emit less radiation than 
they receive, leading to a net warming. In the strato- 
sphere, the lapse rate is negative, leading to net cool- 
ing. Because the height of the tropopause, which sep- 
arates these two layers, is higher at low latitudes due 
to vertical mixing, and decreases abruptly poleward of 
the mid-latitude jet, increasing greenhouse gases will 
enhance the meridional temperature gradient across 
the jet core. This gradient, which is largely zonally 
symmetric, is associated with strengthened westerlies 
within the jet. The strengthened mid-latitude jet al- 
ters the propagation of planetary waves, which can 
potentially feed back on the jet through eddy forcing. 

Any amplification effect would be largest in North- 
ern Hemisphere winter, when planetary wave activity 
is greatest. While we used the November- April cold- 
season for SLP, to facilitate comparison with other 
analyses of the AO, we return now to the conventional 
December-February winter period for easier compari- 
son with model output and with observations of more 
standard meteorological parameters. In that season, 
the enhanced polar night jet strengthens the polar 
vortex over the Arctic [Shindell et at., 1998a; Rind 
et ai, 1998]. Planetary wave refraction is governed 
by wind shear, among other factors, so that enhanced 
wave refraction occurs as the waves approach the area 
of increased wind. In this case, since planetary waves 
are propagating up from the surface, they are re- 
fracted by the increased vertical shear below the en- 
hginced zonal wind. Equatorward refraction of plan- 
etary waves at the lower edge of the wind anomaly 
(Figure 7) leads to wave divergence and hence an ac- 
celeration of the zonal wind in that region. Over time, 
the wind anomaly itself thus propagates downward 
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[Haynes et o/., 1991] from the lower stratosphere to 
the surface. 

The precise location of the enhanced westerlies 
shown in Figure 7 therefore depends upon the inter- 
action between planetary waves, wind shears, eddy 
fluxes and angular momentum fluxes, and is not a 
simple function of the location of the strongest tem- 
perature contrast. This is likely the cause of the dif- 
ference in location between the strongest enhance- 
ment of the temperature gradient, which is at the 
latitude of the jet stream, and the location of the 
strongest zonal wind enhancement, which is at the 
latitude of the polar night jet. Further work will 
be required to fully elucidate the link between the 
meridional structures of the temperature and wind 
responses to external forcing. 

Observations suggest that the Arctic vortex has 
indeed strengthened recently [Tanaka et a/., 1996; 
Zurek et al.^ 1996; Waugh et a/., 1999; Hood et a/., 

1999] . Likewise, a trend towards equatorward wave 
fluxes and downward propagation of wind anomalies 
is also seen in observations [Kodera and Koide^ 1997; 
Kuroda and Kodera, 1999; Ohhashi and Yamazaki, 
1999; Baldwin and Dunkerton, 1999; Hartmann et ai, 

2000] . Other model simulations show qualitatively 
similar effects on planetary wave propagation [Kodera 
et al, 1996; Perlwitz et a/., 2000]. 

The overall equatorward refraction of planetary 
waves and reduced upward propagation at high lati- 
tudes shown in Figure 7 affects the location of wave 
dissipation. The net result is a reduced ability of indi- 
vidual waves to abruptly enhance the residual circu- 
lation and create sudden warmings. In fact, there 
is a strong anti-correlation between the frequency 
of sudden stratospheric warmings and the AO index 
strength [Hartmann et ai, 2000]. 

The behavior of the troposphere is in accord with 
theory stating that angular momentum is in general 
transported in the opposite meridional direction to 
planetary wave energy [Andrews et al, 1987]. In this 
case, angular momentum is then preferentially trans- 
ported poleward, enhancing westerlies, as waves are 
refracted towards the equator in the upper tropo- 
sphere. Geostrophic and hydrostatic balance in the 
atmosphere is maintained by generation of a verti- 
cal circulation cell consistent with the northerly an- 
gular momentum transfer. Increased westerlies are 
transformed into northerly flow by surface friction, 
leading to a cell with rising air in the polar region, 
and descending air at middle latitudes from about 
40°to 55° N throughout the troposphere and the lower 


stratosphere. The effects are clearly visible in the sea- 
level pressure field as a decrease in the Arctic con- 
current with increased mid-latitude pressure. Adia- 
batic expansion of rising air at high latitudes must 
be balanced by radiative heating, which occurs as a 
result of further cooling of the air below its radiative 
equilibrium temperature. Conversely, sinking air at 
mid-latitudes warms. The increased latitudinal tem- 
perature gradient that results is consistent with the 
increased westerly zonal wind around 55°N seen in 
the model and in observations [Baldwin and Dunker- 
ton, 1999; Thompson et ai, 2000]. The surface wind 
anomaly is therefore created by the effect of the in- 
creased lower stratospheric zonal wind on planetary 
waves via linked changes in planetary wave and an- 
gular momentum fluxes. It is this increase in sur- 
face wind that leads to greater advection of warm 
oceanic air over the downstream continents. The role 
of the stratosphere in influencing the behavior of sur- 
face meteorology in the model is in agreement with 
the observed co-variability between the lower strato- 
sphere and the surface [e.g. Perlwitz and Graf, 1995; 
Graf et ai, 1995; Thompson and Wallace, 1998]. It is 
also consistent with the large role played by planetary 
waves in creating the AO pattern itself during mid- 
winter in the Northern Hemisphere, as evidenced by 
the fact that the AO structure amplifies with height 
up into stratosphere only during this season, when 
the zonal flow is conducive to strong wave-mean flow 
interaction [Baldwin and Dunkerton, 1999; Thomp- 
son and Wallace, 2000]. This would also explain why 
there has only been an upward trend in the AO during 
Northern Hemisphere winter. Circulation changes 
arise as a result of the approximately zonally sym- 
metric thermal gradient changes caused by increasing 
greenhouse gases. Planetary wave propagation is a 
secondary response to those changes. Since the ini- 
tial forcing is zonally symmetric, it is perhaps not 
surprising that the response is composed primarily 
of an enhancement of the largely zonally symmetric 
AO pattern. Within this pattern, however, there are 
some longitudinal variations as the centers of wave 
activity are typically downstream of wave generat- 
ing coastlines and mountain ranges. The strongest 
component of the mid-latitude SLP variability is over 
the North Atlantic (the much studied North Atlantic 
Oscillation [e.g. Hurrell, 1995]. Longitudinal varia- 
tions in the model can best be seen by examining 
the driving force for the SLP changes, the middle 
troposphere wave activity flux for stationary plane- 
tary waves [Plumb, 1985]. The model has its max- 
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imum meridional changes occurring over the North 
Atlantic, where there is a distinct equatorward shift 
in wave activity around 500mb (Figure 8), similar to 
the observed change in this sector [Kodera et a/., 1996; 
Ohhashi and Yamazaki, 1999]. This shift has been 
shown to be tied to the AO pattern in other models 
as well [Limpasuvan and Hartmann, 1999]. 

The equatorward movement of the flux conver- 
gence zone allows westerly flow in the northernmost 
Atlantic to proceed more rapidly, while impeding 
flow at more southerly latitudes. Similar correla- 
tions between flow across the North Atlantic and the 
strength of the AO (or NAO) are present in observa- 
tions [Hurrell, 1995; Baldwin and Dunkerton, 1999], 
and indicate a preference for storms to follow the 
more northerly track during the high phase of the 
AO. This shift persists throughout the winter, lead- 
ing to a strong reduction in the energy convergence 
in the region from Iceland to Scandinavia and down 
through the North Sea, which is typically a strong 
convergence zone. Reduced wave energy convergence 
here allows the westerly zonal wind to increase signif- 
icantly, and therefore plays a large role in the down- 
stream warming effects over the Eurasian continent. 
Over East Asia and the Western Pacific, there is 
a westward translation of the wave energy conver- 
gence/divergence pattern, a meridional banded struc- 
ture to begin with. Such a rotation, of course, has lit- 
tle impact on the zonal average zonal flow. There 
is also a small southerly component to the energy 
convergence shift over the Western Pacific. Over- 
all then, there is a tendency toward increased wave 
energy propagation away from the pole, as expected 
from the increased winds induced by the thermal gra- 
dient. This conforms to the generation of an AO type 
response, with SLP decreases in the polar region, and 
increases at middle latitudes, especially over the At- 
lantic and to a lesser extent, the Pacific. It is also 
consistent with a response pattern that is largely zon- 
ally symmetric being created in response to a zonally 
symmetric thermal gradient increase 

5.2 Effect of Arctic ozone depletion 

The depletion of ozone in the Arctic region occurs 
during spring when sunlight returns to the polar re- 
gion. Ozone depletion takes place extremely rapidly, 
so that very soon after depletion has begun temper- 
atures in the lower stratosphere show significant de- 
creases, as there is less ozone to absorb solar ultra- 
violet radiation [Randel and Wu, 1999b]. The cool- 
ing of the polar lower stratosphere from ozone de- 


pletion results in a meridional temperature gradient 
increase, as in the case of increasing greenhouse gases. 
The effects on planetary wave propagation are qual- 
itatively almost identical to those seen in the green- 
house gas simulations (Figure 7). The mechanism 
by which Arctic ozone depletion induces an increased 
AO is thus basically identical to that of the increas- 
ing greenhouse gases, so it is not surprising that the 
two trends both project primarily onto the AO pat- 
tern (Figure 2). However, the Northern Hemisphere 
ozone depletion occurs primarily during late Febru- 
ary, March and April, so that the effect averaged over 
the entire cold-season is much weaker than that in- 
duced by the greenhouse gases (Figure 2), The pri- 
mary effect of severe ozone depletion on the polar vor- 
tex is to make it more persistent into the spring [Graf 
et al, 1998]. 

While there is a relatively steady increase in the 
cold-season averaged strength of the Arctic vortex, 
as shown by the trend in the 30 mb geopotential 
height field, ozone depletion remains quite variable 
during the next four decades since it depends primar- 
ily on the vortex strength only during spring (Figure 
9). Nevertheless, the trend towards a stronger vor- 
tex in the model does lead to significantly enhanced 
ozone depletion relative to what would occur in the 
absence of climate change. This is clearly seen in 
Figure 9, which shows that the maximum predicted 
ozone losses take place well after the maximum abun- 
dances of ozone-destroying halogens are present in the 
stratosphere. 

The ozone recovery trend is fairly constant after 
about 2020 (middle panel), and the increase in the 
AO index also levels off at this time (Figure 3b). The 
vortex, however, continues to strengthen, as shown 
by the geopotential height field (top panel). A closer 
examination of the height index shows a slowdown in 
the rate of increase after about 2020. The saturation 
seen in the surface AO trends, which are transmitted 
downwards by increased planetary wave divergence, 
at around the same time is indicative of the reduced 
influence of planetary wave propagation changes af- 
ter about 2020 in the GCM (note that since there 
are still stratospheric warmings taking place in later 
decades, as demonstrated by the large variability in 
ozone loss shown in Figure 9, it is not obvious why 
the surface trends and planetary wave changes do sat- 
urate). The continued increase in the strength of the 
vortex thereafter is therefore likely driven primarily 
by the radiative impact of the greenhouse gases. The 
continued strengthening of the polar vortex aloft im- 
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plies that the ozone depletion that takes place in the 
mid-21st century will persist later into the spring than 
it would have in the absence of increasing greenhouse 
gases. This effect is overshadowed by the influence 
of the ozone recovery itself on the persistence of the 
vortex [e.g. Zhou et al.^ 2000], through its feedback 
on temperatures, however. This latter effect domi- 
nates the overall persistence, so that the vortex lasts 
longer and longer into the spring as ozone depletion 
increases through about 2020, reaching a maximum of 
about two to three weeks delay in the final warming, 
but then its duration is slowly reduced from about 
2020 to 2040. After this time, it slowly begins to last 
longer again as ozone depletion in the Arctic is mini- 
mal, while greenhouse gases continue to increase 

5.3 Effects of solar cycle variability 

Solar variability impacts the atmosphere initially 
by inducing a thermal gradient in the upper strato- 
sphere [Shindell et a/., 1999b]. In that region, in- 
creased UV radiation at solar maximum relative to 
solar minimum enhances the photochemical produc- 
tion of ozone in sunlit areas. Heating in these al- 
titudes is due primarily to ozone absorption of UV 
radiation. Since there is both more UV and more 
ozone, there is more heating in sunlit areas at solar 
mcLximum than minimum. In the wintertime, the po- 
lar region is in total darkness, so temperatures there 
do not change with solar irradiance variations. Thus 
a latitudinal temperature gradient and a zonal wind 
anomaly are induced in the upper stratosphere (Fig- 
ure 10), and again planetary wave propagation is af- 
fected. The wind anomaly gradually moves down- 
ward [Kodera et al, 1990; Balachandran and Rind, 
1995; Rind and Balachandran, 1995; Shindell et al, 
1999b], and as it reaches the lower stratosphere dur- 
ing the middle part of the winter, the resulting hor- 
izontal and vertical wind shear changes refract plan- 
etary waves equatorward as in the increasing green- 
house gas case. This alters the transport of tropo- 
spheric energy, resulting in temperature changes in 
the lower atmosphere which in turn affect geopoten- 
tial height, causing significant increases at low and 
mid-latitudes and decreased heights near the pole 
[Shindell et al, 1999b]. These height changes are in 
good agreement with an observational record stretch- 
ing back more than forty years, and showing a very 
high level of statistical significance [Labitzke and van 
Loon, 1997], giving credence to the modeled flux and 
temperature changes (a comparison with other mod- 
els is given in Shindell et al [1999b]). The zonal 


mean zonal wind and tropospheric EP flux changes 
in January in fact look fairly similar to those seen in 
the increasing greenhouse gas simulations. December 
and February do not, however (Figure 10), as waves 
fluxes in the increasing greenhouse gas simulations 
show little month-to-month change in the upper tro- 
posphere/lower stratosphere region. 

In contrast to the case of increasing greenhouse 
gases, the initial temperature changes induced by so- 
lar variability take place much higher in the strato- 
sphere. The troposphere sees no direct change in 
thermal gradients, but instead sees altered planetary 
wave propagation as the driving force behind changes. 
It is the changes in wave propagation that lead to 
the altered temperatures and circulation patterns in 
the lower stratosphere, as well as in the troposphere. 
Planetary waves are generated by topography and 
land-ocean contrasts, so are by nature longitudinally 
asymmetric. Thus the impact of solar forcing in the 
lower stratosphere is not as symmetric as that of the 
greenhouse gases, though the zonal mean changes are 
at times broadly similar. Stratospheric water vapor 
also causes large changes in the upper stratosphere 
[Shindell, 2000], and the wave fluxes in the run with 
increasing water appear to be a combination of those 
seen in the solar cycle runs and those seen in the 
greenhouse gas simulations, though primarily the lat- 
ter. 

At solar maximum, averaged winter wave activity 
flux at 584 mb again shows an overall equatorward 
shift relative to solar minimum, but in this case, the 
main activity takes place over the Western Pacific. 
Over the North Atlantic, there is actually a slight 
poleward shift in overall wintertime wave energy con- 
vergence zones, but the pattern varies strongly over 
the course of the winter. There is little signal in De- 
cember, while the February pattern looks the most 
similar to the greenhouse gas induced pattern, with 
a slight southerly shift in the North Atlantic energy 
convergence. Over the Western Pacific, there is a 
distinct eastward and equatorward shift in wave en- 
ergy convergence, which is stable during the entire 
winter and larger than the shift seen in the green- 
house gas runs. This large, more southerly energy 
propagation is responsible for an overall zonal mean 
equatorward shift in response to solar cycle forcing 
[Shindell et al, 1999b], just as in the greenhouse gas 
simulations' zonal mean responses. This greater wave 
energy divergence at high latitudes allows the zonal 
wind to accelerate in that region, creating a strong in- 
crease in the westerly flow across the northern-most 
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Pacific. These results are consistent with a response 
that begins with a zonaily symmetric thermal gra- 
dient increase in the upper stratosphere, but that 
takes place primarily via altered wave propagation 
once it reaches the lower atmosphere. The centers of 
activity are therefore downstream from the areas of 
largest wave activity, i.e. the eastern coasts of Asia 
and North America (neglecting the orographic forc- 
ing from the Himalayas at lower latitudes). The wave 
energy flux changes affect the downstream synoptic 
pressure systems as well as the overall westerly flow, 
so that the response pattern is much less zonaily sym- 
metric than with the greenhouse gases. The pattern 
therefore projects onto many EOFs of sea- level pres- 
sure (Figure 2). Note, however, that the use of fixed 
sea-surface temperatures in these runs may have lim- 
ited the model’s response at the surface. 

The increased surface winds in the solar experi- 
ment reflect a strengthening of the Pacific subtropi- 
cal and Azores high pressure systems, and so the in- 
creased westerlies over both ocean basins are merely 
the northward portion of an increased anticyclonic ro- 
tation. The impacts on downstream surface warming 
are significant, as shown in Figure 11. Solar forcing 
leads to downstream warming at higher latitudes than 
do greenhouse gases. Differences are especially stark 
over Asia, where the greenhouse gases lead to a warm- 
ing over most of the continent north of 45 degrees, 
but in the solar experiments, the return flow down 
from the Arctic leads to cooling south of 60 degrees. 
The influence of using fixed sea-surface temperatures 
is large, however, so further work in this area will be 
necessary to obtain more reliable results 

5.4 Effect of Volcanic Eruptions 

Volcanic forcing induces SLP variations that occur 
predominantly via the AO pattern, though not as ex- 
clusively as in the greenhouse gas experiments or as 
in the past three decades of observations (Figure 2). 
The enhanced AO leads to a wintertime continental 
warming following eruptions, as clearly seen in ob- 
servations [e.g. Robock and Mao^ 1995; Kelly et o/., 
1996]. Changes in wave propagation are quite similar 
to those seen in the greenhouse gas case, especially in 
the mid to high latitude troposphere (Figure 12). The 
injection of stratospheric aerosols causes enhanced ab- 
sorption of shortwave radiation, leading to heating in 
the sunlit lower stratosphere [e.g. Rind et oi., 1992; 
Graf et al, 1993]. Along with aerosol- induced long- 
wave cooling at high latitudes, the volcanic forcing 
thus creates an increased meridional temperature gra- 


dient in the lower stratosphere, particularly in winter. 
This induces a positive zonal wind anomaly, which 
then leads to equatorward refraction of upward prop- 
agating planetary waves. Similar behavior is seen in 
observations [Kodera, 1995]. Note that the magni- 
tude of the AO index response found here, -0.4 mb 
between volcanic and non-volcanic years, is based on 
a fairly small average aerosol loading. As noted, the 
aerosol loading immediately following the eruptions 
of El Chichon and Mt. Pinatubo were 4 and 6 times 
larger, respectively, than the ’volcanic years’ average 
used here. The model’s AO index increase after these 
eruptions was much smaller than what a linear ex- 
trapolation would yield, implying that saturation ef- 
fects take place in this forcing as well as in the green- 
house gas forcing. 

The importance of the stratosphere to the response 
seems to be fairly similar in other models. A large 
set of tropospheric models used in the Atmospheric 
Model Intercomparison Project (AMIP) could not re- 
produce the observed wintertime warming over the 
Northern Hemisphere continents after El Chichon ’s 
eruption, which can most likely be attributed to de- 
ficiencies in those model’s simulation of stratospheric 
processes [Mao and Robock, 1998]. In contrast, mod- 
els with a somewhat more complete representation 
of the stratosphere do show an enhanced AO type 
pattern, as in the GISS stratospheric model, with 
significant continental surface air temperature in- 
creases during wintertime following eruptions [Kirch- 
ner et ai, 1999; Robock et al, 1999]. It is interest- 
ing that a simulation of the response to the erup- 
tion of Mt. Pinatubo, using a version of the ECHAM 
model with a top at 10 mb, does seem to reproduce 
the observed magnitude of AO enhancement [Kirch- 
ner et a/., 1999]. (Note that though the model top is 
in the middle stratosphere at 10 mb, this model nev- 
ertheless has a much more realistic stratosphere than 
the AMIP models). This is likely related to the ex- 
tremely large forcing of the Pinatubo aerosols, which 
warmed the lower stratosphere by roughly 3 degrees. 
This induces a meridional thermal gradient as large 
as that induced by increasing greenhouse gases over 
about 50 years. Since the greenhouse gas forcing sat- 
urates in the GISS model after increasing for roughly 
60 years [Shindell et al, 1999a], presumably as the 
polar vortex has become strong and stable enough so 
that nearly all planetary waves are already refracted 
away, any further forcing is not effective. So while a 
model with a top at 10 mb might not fully capture 
alterations in planetary wave propagation, amplifica- 
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tion of such a large, initial thermal gradient may not 
be necessary to produce the appropriate response in 
this particular case 

6. Conclusions 

A series of general circulation model simulations 
demonstrate that changes in the stratosphere can 
influence surface climate. Either anthropogenic or 
natural perturbations to the stratosphere are capa- 
ble of causing changes to the stratospheric thermal 
structure, which in turn affect stratospheric circula- 
tion. These wind changes alter the wind shear, affect- 
ing propagation of planetary waves out of the tropo- 
sphere so as to amplify the lower stratospheric wind 
changes during Northern Hemisphere winter. Fur- 
thermore, the deflection of upward propagating waves 
takes place at steadily lower and lower altitudes, and 
so exerts a significant impact on surface climate. This 
fundamental role for planetary waves, which are rel- 
atively absent in the Southern Hemisphere, is consis- 
tent with studies showing that the radiative impact 
of ozone depletion is able to account almost com- 
pletely for observed wintertime temperature trends 
in the Antarctic, but not in the Arctic [Randel and 
Wu, 1999b]. 

The model simulations have shown that a surface 
climate response is induced by increasing greenhouse 
gases, ozone depletion, volcanic eruptions, and solar 
variability. Except for solar variability, the others all 
affect surface climate at Northern Hemisphere middle 
and high latitudes during winter primarily by favoring 
the positive phase of the AO, the dominant natural 
variability pattern, while concurrently strengthening 
the polar vortex aloft. Nonetheless, only increasing 
greenhouse gases seem capable of causing an AO trend 
in our model as large as the observed trend. 

Other studies [Fyfe et aL, 1999; Paeth et al., 1999; 
Zorita and Gonzalez- Rouco^ 2000] note that some 
models lacking detailed resolution of stratospheric 
dynamics, exhibit non-zero AO/NAO trends when 
forced with increasing greenhouse gases. Similar mod- 
els, such as the ECHAM3 [Perlwitz et al, 2000] and 
HadCM3 [R. MacDonald, personal communication, 
2000] also exhibit non-zero trends that are nonethe- 
less weak compared to the magnitude observed in re- 
cent decades. This is consistent with the behavior 
of the GISS model described above, where only the 
inclusion of stratospheric dynamics results in an AO 
trend of magnitude comparable to that suggested by 
observations in recent decades. 


The AO sensitivity of models developed by several 
groups, including GISS, which lack realistic represen- 
tations of the stratosphere, is much weaker than the 
GISS stratospheric model. It seems that models with 
an overly simplistic stratosphere fail to adequately 
capture the positive feedback of planetary wave flux 
changes on circulation anomalies seen in the strato- 
spheric model. Furthermore, the simulated trends in 
those models are not primarily composed of the AO 
pattern (with the exception of the ECHAM4/OPYC 
model, A. Robertson, personal communication, 2000), 
as are the observations, while the stratospheric model 
does reproduce this behavior. 

We hope that future studies will present quantita- 
tive trends of the AO index (along with its statisti- 
cal significance) for comparison to the observed values 
and other modeling results (as opposed to the normal- 
ized 'AO/NAO’ indices), while quantifying the extent 
to which their trend in SLP resembles the AO spatial 
pattern. In the meantime, the behavior of the tropo- 
spheric models cited above remain in contrast to the 
GISS stratospheric model. As yet, only the latter has 
shown an AO trend of roughly the right magnitude 
and with the same predominance in the SLP trend 
as the observations. It remains to be demonstrated 
that a model lacking stratospheric dynamics is able 
to match both these aspects of the observations. 

The surface climate response of the GISS strato- 
spheric model to the injection of volcanic aerosols 
into the stratosphere corresponds well with obser- 
vations. This provides important evidence that the 
modeled linkage between forcings and the AO does 
indeed occur in the atmosphere. While longer-term 
AO changes will need to continue for many more years 
to be fully accepted as a distinct trend, the AO en- 
hancement after eruptions is already quite clear in 
the observational record. Since the AO responses to 
greenhouse gas increases, ozone depletion, and vol- 
canic aerosols all take place via the same mechanism, 
the clearly observed response to large volcanic erup- 
tions suggests that the same linkage is likely at work 
in the case of the slower, longer-term greenhouse gas 
and ozone depletion forcings as well. 

These results support our earlier conclusion that 
much of the trend in the AO and in the Arctic vor- 
tex strength over the past few decades should be at- 
tributed to increasing greenhouse gases, and that the 
trends are likely to continue. 
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Figure 1. Surfeice temperature trends (C). Top left: Observed surface temperature trend from 1965 to 1995, 
December to February average [Hansen et a/., 1999]. Top right: The AO contribution to the total trend, obtained 
by regressing the temperature trend onto the AO spatial pattern as in Thompson et al. [2000]. Bottom center: 
The AO component of the total trend in the GISS stratospheric model, similarly obtained by regressing the 
model’s temperature trend onto its AO spatial pattern. Grey areas are those where no data is present. Note 
that the difference in grey areas between the top left and top right panels results from the use of slightly different 
temperature data sets, which do not match exactly in areas with few observations (Arctic Ocean and low latitude 
Atlantic and Pacific). 
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Percent of Northern Hemisphere Winter SLP C h a ng e in each EOF 




Figure 2. Percentage of the total Northern Hemisphere November-April sea-level pressure change in each EOF. 
Trends from e^w;h indicated model simulation were projected upon the control run EOFs. The magnitude of the 
trend or change in the AO index induced by each forcing is also given. Values are calculated using a linear least- 
squares regression through the observations or model output. Uncertainties are given in Table 1, and are 95% 
confidence limits based on a student’s-T test, assuming each winter is independent (Note that assuming red noise 
instead of white has a minimal effect on these results). Observations are updated from TYenberth and Paolino 
[1980]. 
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Figure 3. Northern Hemisphere sea level pressure trends. The left column shows the leading EOF (the AO), 
while the right column shows the principle component (PC) for: the observations (top), the stratospheric model 
(middle), and the ’tropospheric’ model (bottom). The PCs are scaled to have units of mbar and the spatial pattern 
has unit amplitude poleward of 60® N. The red curve is a 10-year smoothed fit to the PC. Note that this figure is 
similar to Figure 2 of Shindell et al [1999a], but that in that figure the PC of observed SLP was inconsistently 
scaled with respect to the model PCs). 
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Zonal wind change over the North Atlantic 
(48N-64N.4SW-5E) 






Figure 4. Observed and simulated zonal wind change over the North Atlantic (48°N - 64*^N, 45^W - 5°E). Positive 
values are westerly winds. Values are computed as in Figure 2, except that the entire period of the NCEP reanalysis 
is used (1951-1997), compared with the trend over the entire length of the model simulations. 

Surface Air Temperature Change over Sberia 
(48 N - 56 N, 62 E - 80 E) 





Figure 5. Observed and simulated surface air temperature change over Siberia (48°N 56°N, 62°E 80°E). 

Vcdues are difference between the average over this location and the 48°- 56°N zonal mean in order to separate the 
regional signal from the latitudinal average. Values are computed as in Figure 4. 






Surface Air Temperature Change over Newfoundland 
(48N-56N.60W-80W) 
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Figure 6. Observed and simulated surface air temperature change over Newfoundland (48°N 56°N, 60°W - 

80°W). As in Figure 4, values are difference between the average over this location and the 48°- 56°N zonal mean 
in order to separate the regional signal from the latitudinal average. Values are computed as in Figure 4. 
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Greenhouse Gases (DJF) NCEP Analysis (DJF) 





Figure 7. Trends in zonal mean zonal wind (m/s) and planetary wave activity (m^/s^) in the greenhouse gas 
simulations (left) and in analyzed observations. Model trends are shown for the five decades during which the AO 
increased steadily between the initial spinup and eventual saturation. Observed trends are differences between the 
high AO years 1981, 1983, 1984 and the low AO years 1979, 1982, and 1992 using NCEP analyses of observational 
data. Values are the sum of the total trend components in the first two EOFs of these trends, as given in Ohhashi 
and Yamazaki [1999], figures 7 and 8 (in a more recent paper, Hartmann et al. [2000] show similar results). Since 
the first two EOFs contain nearly all the variability, these should be comparable to the model’s total trend values. 
Vertical fluxes have been scaled by the inverse of pressure for convenience in all plots. The maximum vector length 
corresponds to 7 m‘^/s^. The total AO index change in the model over this five decade period was about 4 mb, 
somewhat larger than the 1979 to 1992 variability of about 2.5 mb. 

Horizontal Rux Changes. GHGs and Obs (m/s*) 



Figure 8. Horizontal wave ctctivity flux trends (m^/s*^) at 584 mb in the model and at 500 mb in observations. 
Model results are as in Figure 6. Observations are the difference between the five year means of 1989-1993 and 
1965-1969, providing an estimate of the long-term trend (identical to Fig. 12b, Kodera et al [1996]). The maximum 
vector length corresponds to 24 m‘‘^/s*^. 
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Year 

Figure 9. TVends in the strength of the Arctic vortex, Arctic ozone depletion, and stratospheric chlorine loading. 
The top panel shows the cold-season (November- April) averaged principle component of the leading EOF of 30 mb 
geopotential height from the stratospheric model with increasing greenhouse gases only, with the 10-year running 
mean displayed as a thicker line. The middle panel shows minimum column ozone (DU) north of 65°, averaged over 
the last three days of March, as seen in TOMS data through 2000 (filled circles) and in the model (open circles). 
Note that model vsJues do not include intrusions of low-ozone air from lower latitudes into the Arctic, so they tend 
to be on the high side of observed values. The bottom panel shows the past and projected lower stratospheric 
chlorine lo 2 wling used in the model (independent of height or latitude) based on World Meteorological Organization 
[1999]. TOMS data obtained from NASA Goddard Space Flight Center via http : // jsocky . gsf c . nasa . gov. While 
ozone loss is often greatest at the end of March, in 2000, ozone column values below 250 DU were seen in early 
March, but did not last until the end of the month. 







Figure 11. Latitude versus longitude December- February surface temperature changes between solar maximum 
and minimum (top), and from a run with increasing greenhouse gas^ (bottom). Small values are white. Note that 
the solar cycle simulations used fixed sea-surface temperatures, limiting their response. 
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Figure 12. Differences in zonal mean zonal wind (m/s) and planetary wave activity (m‘'^/s‘'^) between volcanic and 
background years. Maximum vector length corresponds to 3.4 m /s . 










